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ABSTRACT 

We critically analyze the Bell et al. findings on "anomalous" widths of high-order Hydrogen radio 
recombination lines in the Orion Nebula at 6 GHz. We review their method of modified frequency 
switching and show that the way this method is used for large An is not optimal and can lead to mis- 
interpretation of measured spectral line parameters. Using a model of the Orion Nebula, conventional 
broadening theory, and Monte Carlo simulation, we determine a transition zone n = 224, . . . , 241 
{An = 11, . . . , 14), where measurement errors grow quickly with n and become comparable with the 
measurement values themselves. When system noise and spectrum channelization are accounted for, 
our simulation predicts "processed" line narrowing in the transition zone similar to that reported by 
Bell ct al. We find good agreement between our simulation results and their findings, both in line 
temperatures and widths. We conclude, therefore, that Bell et al.'s findings do not indicate a need to 
revise Stark broadening theory. 

Subject headings: H ii regions - line: profiles - methods: data analysis - methods: statistical - radio 
lines: ISM - techniques: spectroscopic 



1. INTRODUCTION 

In the 1990s, Morley Bell and coauthors developed a 
technique for measuring weak spect ral lines by reducing 
broad baseline variations ()Bellll 19971 ) . This technique was 
referred to as "modified frequency switching" and was 
used to detect of weak atomic and molecular /maser lines 
in spectra of gaseous nebulae, circums tellar envelopes, 
and star formation regions (see, e.g., I Bell fc FeldmanI 
flMllBelTet al.lll992t [l993^. 

Bell and coauthors used this technique to measure 
radio recombination lin e (RRL) widths and tempera- 
tures (|BeU et all [2000l hereafter BASV). At 6 GHz, 
they found that the "...[processed H] lines at large n 
are both narrower and stronger than expected from the- 
ory..." and suggested, ^^This behavior is ... inconsis- 
tent with Griem's theory ...". This publication, with 
its subtitle: "Confrontation with theory at high princi- 
pal quantum numbers" induced a wave of publications 
where BASV's finding was called an " anomaly" , "puz- 
zle" and even "mystery". For example, lOkj [2004l paper 
is titled "On the puzzle o f the ob served narrowing of ra- 
dio recombination lines," lGriera ,2005 concludes his pa- 
per by writing that the result "...remains a mystery...," 
iGavrilenko fc Oksl 12007 calls the result a "dramatic dis- 
crepancv." and lWatsonI 2006 concludes his abstract with 
"Thus this mystery is not resolved by the present calcu- 
lations" . 

While some authors have sought an explanation of 
these find ings in th e revision of S tark broadening the- 
ory (e.g., I0ksll2003 . lWatsonll2006[ ). some remain skep- 
tical about the modified frequency switching technique 
or suggest that it require s verification. For example, 
Ivon Prochazka et al.ll2010l writes "...it is easy to distort 
the RRL line sh ape using dat a reduction techniques (Bell 
et al. 2000)." iGriemI 120051 referred to private commu- 
nication with Bell stressing that "...measurement errors 
related to large reduction factors from unprocessed to 
processed Voigt profile widths can probably not be ex- 
cluded (M. B. Bell 2004, private comm unication)". In 
their book. lGordon fc Sorochenkoll2009l concluded: "The 



Bell et al. (2000) results are so different from what had 
been expected, and the observing technique is so new, 
that we suggest waiting for an independent confirmation 
of the observations before accepting a fault in the present 
theory of RRL Stark broadening." 

BASV's findings are based on frequency switched or 
"process ed" observationa l data, which were re cently pre- 
sented in lBell et al.ll20lil and re- interpreted in lBellll201lL 
Their findings result from recursive frequency switch- 
ing, in software, of observational data that were ini- 
tially recorded at the telescope using hardware frequency 
switching. 

In this paper, we investigate the frequency switching 
technique and demonstrate that, if applied correctly, it 
has a number of advantages. Frequency switching re- 
moves many gain variations and does not require subjec- 
tive estimates of the zero level of a spectrum. Though 
the technique helped BASV to detect RRLs with An 
greater than An — 6 (Smir nov et al. 1984 ). the way the 
method was used was not optimal and, as such, cannot 
be used to test the theory of spect ral line broadening. 
We p resent simulations based on the lLockman fc Brownl 
119751 model of the Orion Nebula and conventional the- 
ory of spectral line broadenin g. We apply obs ervational 
specifications from BASV and lBell et al.l[20TTl including 
frequency range, channel width, frequency switching off- 
set, number of frequency switching overlaps, and noise 
temperature rms. Results of our simulation demonstrate 
good agreement with BASV's findings, both in line width 
and temperature. The computed "processed" widths ex- 
hibit narrowing similar to that reported in BASV. We 
show that BASV's spectral line "narrowing" is the re- 
sult of the way the observational data were processed 
and that BASV's findings do not contradict the existing 
Stark broadening theory. 

2. SOFTWARE FREQUENCY SWITCHING 

In this section, we introduce the data reduction tech- 
nique used by BASV to systematically acquire informa- 
tion about weak s pectral fea tures in the presence of base- 
line fiuctuations. iBelll 119971 established and named this 
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technique SOFMOR ("small-offset frequency-switching 
multiple overlap reduction"). This technique was then 
called "modified frequency switching" (MFS) in BASV. 
SOFMOR and MFS are frequency switching techniques 
that use a freq uency o ffset fixed in the receiver hard- 
ware (see, e.g., iRobins on 1964), followed by recursive 
frequency switching in software. Here, we refer to this 
technique (whether conducted in hardware or software) 
as FS (frequency switching). 

In FS, the original spectrum is subtracted from a copy 
which is shifted by a number of channels (an offset). 
Mathematically, it can be related to a finite difference 
iteration. 

If the offset is h, then the notation 

fp^fix+ph) (1) 

can be used. Formulae for the results of m FS iterations 
(m finite differences) are then 
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and so on (see, e.g., iBeveJUOSTL p. 449: IZwillingeJ[2002L 
p. 705), therefore 
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When /(x) is given with uncertainty ±ao, where ctq 
is standard deviation, the resulting uncertainty in 
can be deriv ed from the error propagation rule (see, e.g., 
lTavlori[l997l p. 75) as 
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Given that an FS spectrum overlaps with itself, the case 
of dependent samples must be applied (see below). 

If f{x) is a single channel feature (digital analogue of 
(5- function), that is 



f{x) = S{xo) 



1, X = Xq 
0, otherwise, 



following Equation ([2]), m FS-overlaps generates a se- 
ries of equidistant peaks with amplitudes (™)(— 1)™^'^- 
This is illustrated in the following Pascal's triangle for 
six consecutive overlaps or three cycles: 
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Each cycle of two overlaps results in a symmetric pat- 
tern with a prominent central feature. The amplitude of 
the central feature after i FS-cycles is the central bino- 
mial coefficient: 



I A: I = 



= 2, 6, 20, 



{i = l, 2, 3,...) . (4) 



According to Equation Q, the standard deviation af- 
ter i FS-cycles 




= \/6, VtO, 7924, 



resulting in a signal-to- noise ratio (S/N) of 

/'2i^ 

= 0.82, 0.72, 0.66, ... 



(i = l, 2, 3, 4,. 

(5) 




{i^l, 2, 3,.. 



Using Stirhng's formula (see, e.g.. lAbramowitz fc Stegunl 
119651 ) we obtain 



1/4 



S/N, ~ - = 0.89 



(7) 



Therefore, BASV's S/N tends to decrease with the num- 
ber of FS-cycles, i, as . 

When independent samples are superimposed, the 
standard deviation is scaled as ct; = 2, 4, 8, . . . , which 
results in the growth of the S/N: S/N, = Ci) f^' = 
1, 1.5, 2.5, ... (i = 1, 2, 3, ...). This common 
approach is not applicable to the FS technique used in 
BASV. Instead, they overlap the once recorded spectrum 
with itself and therefore the case of dependent samples 
must be used, as derived in Equations ©-([T]). 

Though the first overlap in BASV occurs in hardware, 
where two independent samples are mixed, the S/N re- 
mains unchanged because the first overlap of two samples 
produces a \/2 increase in the resulting standard devia- 
tion regardless of whether the two samples are dependent 
(software FS) or independent (hardware FS). 

Figure [T] shows a Gaussian feature characterized by its 
full width at half-maximum (FWHM), w, and its trans- 
formation after 2, 4, and 6 FS-overlaps are applied to 
the initial feature. The frequency offset or FS-offset, h, 
is chosen to be much greater than the width, h ^ w. 
Additional features appear offset from the central fea- 
ture, whose offsets are multiples of h. All the features 
are shifted relative to the position of the original feature, 
Xq. After i FS-cycles, the central feature is located at 
X = Xq — i h. These shifts are removed in Figures [T] and 
m such that all features are centered vertically. Addi- 
tionally, all features are normalized and inverted (where 
necessary) to compare their shapes and widths. 

li h ^ w, the additional features (referred to as "ref- 
erence images" in iBell 1997} can be removed by adding 
or subtracting scaled versions of the central feature to 
them — a process referred to as "cleaning" in iBelll [T9971 
For example, after two FS-overlaps, the additional fea- 
tures shown in Figure [1] can be removed by adding the 
central feature divided by 2. Dividing this result by —2 
recovers the original feature exactly. After four overlaps, 
shown in Figure [1] the additional features are removed 
by adding —1/6 and 4/6 of the central feature to them. 
Dividing the central feature by 6 recovers the original 
feature exactly. 

This normalization procedure keeps the central fea- 
ture's amplitude unchanged regardless of the number of 
FS-cycles, i. However, according to Equation ([7]), the 
noise level grows as i^^'^ . 

NRAO's 140 ft telescope was used by BASV to observe 
the Orion Nebula at 6 GHz. Based on the 140 ft data 
archive of the 1992 April observation (integration time of 
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Figure 1. Transformation of a Gaussian sliape after 2, 4, and 6 
FS-ovcrlaps. Tlie initial width, w, is much less than the FS offset, 
h; w/h = 0.2. This plot is a graphical depiction of Pascal's triangle 
in Section [2] To compare the shapes and widths, wc invert (where 
necessary) the resulting spectra and normalize the peak intensity 
of the central feature to the intensity of the initial Gaussian. 




-I 1 , , , — I — , , , 1 

012345678 
Unprocessed width / offset 




A 2 overlaps 




4 overlaps 



6 overlaps 



Frequency 

Figure 2. Transformation of a Lorentzian shape after 2, 4, and 6 
FS overlaps. The initial width, w, is greater than the offset, h. The 
ratio = 4 is used, which corresponds to BASV's An = 20 case. 
To compare the shapes and widths, we invert (where necessary) the 
resulting spectrum and normalize the peak intensity of the central 
feature to the intensity of the initial Lorentzian. 

~48 hours, system temperature of ~120K, and channel 
width of ~ 78 kHz) , a noise temperature Trms ~ 1 mK 
could be achieved. Given that BASV used six overlaps 
(one hardware and five software FS-overlaps) or three cy- 
cles, Equation ([7]) implies an "frequency-switched" tem- 
perature rms of Trms ~ 1-5 mK. 

Figure [2] illustrates the case when the FS-ofFset is less 
than the original line width, h < w, specifically h — w/A. 
This figure shows the transformation of a Lorentzian fea- 
ture after 2, 4, and 6 FS-overlaps. To compare the shapes 
and widths, we normalize the peak intensity of the cen- 



Figure 3. Reduction factors (processed over unprocessed) for 
the peak line temperature (top plot) and the line width (bot- 
tom plot) as a function of the normalized unprocessed width, 
w/h. Both plots are calculated for six FS-overlaps. Vertical lines 
show minimal and maximal values of w/h for BASV's observations: 
f ^ w/h < 4.5. Crosses show reduction factors predic ted by the 
three-component Orion Nebula model of lLockman & Browni/19751 . 
Horizontal dashed lines in the bottom plot are the asymptotes for 
Gaussian and Lorentzian curves as w/h — ^ oo. These plots demon- 
strate that, for a fixed FS-offset h, significant reduction in spectral 
line intensity and width results from FS as the unprocessed line 
width increases due to Stark broadening. 

tral feature to the intensity of the initial feature. Figure 
[2] shows that the central feature of the processed line is 
much narrower than the width of the original Lorentzian. 
Given that FS processing reduces line widths and ampli- 
tudes when h < w, BASV refer to them as "processed" 
widths and amplitudes. 

The Voigt profile, which models the RRL shape, 
is a convolution of Gaussian and Lorentzian profiles, 
where the Gaussian represents the Doppler-broadened 
contribution and the Lor entzian represents the i mpact- 
broadened contribution (jGordon &: Sorochenkol [2009.) . 
Figure [3] shows how the amplitudes and widths of FS- 
processed Gaussian and Lorentzian profiles vary as a 
function of the normalized unprocessed width, w/h. 
These plots show amplitude and width reduction factors 
(processed to unprocessed ratio) calculated for six FS- 
overlaps. The vertical lines show the limits of BASV's ex- 
periment: 1 < w /h S, 4.5. Crosses a.re red uction factors 
predicted by the iLockman fc Brownlll975l Orion Nebula 
model (see the next section). Horizontal dashed lines in 
the bottom plot are the asymptotes for Gaussian and 
Lorentzian cases as w/h — >■ oo. 

3. APPLICATION OF FS TO MODEL SPECTRA 

The frequency switching technique used in BASV con- 
sists of two independent steps: (1) overlapping the spec- 
trum and (2) "cleaning" the overlapped spectrum. Over- 
lapping six times creates a set of "reference images" for 
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Figure 4. Processed width vs. n for different number of FS over- 
la ps at 6 GHz. Sim ulation is conducted for the Orion Nebula model 
of lLockm a n fc Brown , 1975: no noise is taken into account. Dashed 
line shows BASV's FS-offset (eight channels = 31.2 km s'^). The 
upper (thick) curve shows the unprocessed width w vs. n. The 
lower curve corresponds to the BASV's case of six FS-overlaps. 
This plot demonstrates that when h < w, the processed width is 
increasingly insensitive to the significantly changing unprocessed 
width as the number of overlaps increase. This indicates that the 
FS technique in BASV's case is a poor probe of line broadening. 

each line feature in the spectrum. The result of over- 
lapping (step 1) is illustrated in Figures [1] {h ^ w) and 
[2] (/j < w). When h < w (BASV's case), overlapping 
reduces the line widths to below the FS-ofFset, h. 

After overlapping, the reference images are removed. 
In Bell (1997) and BASV, this "cleaning" procedure was 
done manually, one line at a time, starting with the 
strongest line and working to the weakest. BASV's re- 
ported line widths after cleaning are in the range of 18 to 
30 km s^^, which is less than or about the Doppler width 
and less than the FS-offset (Figures 2 and 4 in BASV). 
The conclusion made in BASV about line "narrowing" is 
based on interpretation of these "processed" widths. 

A multi-component non-LTE radiative transfer simu- 
lator was created to explore the effect of FS on model 
RRL spectra. For the Orion Nebula we use the three- 
component model o flLoc kman & Brown 1975, which sim- 
ulates the beam size of the 140 ft radio telescope. The 
model consists of a compact (0.043 pc) dense (A^e = 
lO'''^ cm~'^) symmetric core located behind two extended 
(0.56 and 2.50 pc) layers of gas of lower density {Ne — 
10^'^ and lO^ "^ cm~'^). The corresponding temperatures 
are Te =7500 K (core), 10,000 K, and 12,500 K (outer 
layer) . This model predicts continuum emission and Una 
and Hn/3 spectr a in agreement with obser vations between 
1 and 90 GHz (iLockman fc Brownl ll975f). 

Figures HHHl show the model results for processed 
line temperatures and processed widths versus principal 
quantum number, n, for the Orion Nebula model at 6 
GHz. To avoid ambiguity associated with the cleaning 
procedure, we process each modeled spectral line sepa- 
rately. The processed line width is the FWHM of the 
central feature measured above the zero level (see Figure 

Figure S] shows the processed widths (FWHM) for 2, 
4, and 6 FS-overlaps (1, 2, and 3 cycles). Following 
BASV, we used six overlaps and an offset of 625 kHz = 
31.2 km s~^, which corresponds to eight spectral chan- 
nels (channel width of 78.13 kHz = 3.9 km s~^). Noise 
was not added to the model spectra when computing 
these curves. The dashed line shows BASV's FS-offset 



Figure 5. Line temperature, T^, vs. n for the Orion Nebula at 6 
GHz. Open circles represent line temperatures obtained by BASV: 
processed line temperatures after six FS-overlaps (Bell et al. 2011, ). 
The theoretical curves result from Lockman & Brown 1975 model. 
Upper and lower solid curves show unprocessed and processed line 
temperatures from the model spectra without noise. Upper and 
lower dashed curves show processed line temperatures from model 
spectra with 3 and 1 mK noise, correspondingly. Vertical lines are 
positioned at n values corresponding to theoretical processed tem- 
peratures = 3mK (left) and Ti^ = ImK (the right vertical line). 
This plot indicates that our model of processed line temperatures 
agrees with BASV's findings for (n. An) < (224, 11). Above this 
limit, our model suggests BASV's results are dominated by noise 
fluctuations. 




Figure 6. Processed width (six FS-overlaps) vs. n computed for 
different FS-offsets (8, 14, and 20 channels) for the Orion Neb- 
ula model at 6 GHz. Dashed curves — no noise; solid curves — 1 
mK noise added to model spectra. Open and filled circles show 
the processed widths obtained by BASV for the Orion Nebula and 
W51. The upper thick solid curve shows the unprocessed width 
vs. n. Vertical lines are positioned at n values corresponding to 
theoretical processed temperatures = 3 mK and Tj^ = 1 mK 
(see Figure [Sjl. The lower thick solid curve corresponds to BASV's 
eight-channel FS-offset and 1 mK noise fluctuations. It agrees 
well with BASV's processed findings for (n. An) < (241,14). In 
the transition zone (between the two vertical lines), the modeled 
width behavior imitates "narrowing" of the processed spectral lines 
reported by BASV. 

(eight channels — 31.2 km s~^). After six FS-overlaps, 
the processed width (FWHM) remains at the level of 
^70%, . . . , 85% of the offset, reaching a maximum at 
n ~ 260, then drops down to 81% at higher n. This be- 
havior reflects a growing contribution to the line shape 
from the nebula's low-density outer layer at very high n. 
The greater contribution of the outer low-density layer at 
high n also explains the deviation of the reduction factors 
(crosses) in Figure |3] from the "Lorentzian" behavior. 
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Figure [5] shows the hue temperature, Tl, versus n 
for the Orion Nebula. Open circles are BASV's pro- 
cessed Tl. Solid curves show unprocessed (upper curve) 
and processed for the model spectra without noise. 
Dashed curves show processed T]^ for model spectra with 
1 and 3 mK noise added. These curves are obtained us- 
ing the following Monte Carlo method. For each modeled 
spectral line feature, normally distributed noise is added, 
then FS-overlapped six times. Tl and w of this processed 
feature are then estimated by a Gaussian fit, as done by 
BASV. The least-squares-fitting algorithm searches for 
the best fit within a ±2 channel interval about the cen- 
tral line frequency. These steps were trialed 1000 times 
to determine the mean and standard deviation values of 
Tl and w for each processed spectral line feature. This 
method was carried out for the expected 1 mK noise rms 
and, additionally, a 3 mK noise rms. Two vertical lines 
in Figure [5] are positioned at n values corresponding to 
the processed, noise- free model temperatures Tl — 3 mK 
(left line) and Tl = 1 mK (right line). Given BASV's 
Trms ~ 1 mK, the vertical lines represent the 3a and la 
limits. The corresponding limiting values of n = 224 
(An = 11) and n — 241 (An = 14) are also indicated by 
vertical lines in Figure [6l 

Processed widths versus n for 8, 14, and 20 channel off- 
sets are shown in Figure [S] for the 6 GHz Orion Nebula 
model. For comparison, the thick upper solid curve is un- 
processed noise-free model widths. Dashed curves in Fig- 
ure [6] correspond to processed noise-free model widths, 
while solid curves are the corresponding processed widths 
resulting from the Monte Carlo method described above 
using 1 mK noise added to the model spectra. Open and 
filled circles in Figure |6] are BASV's pr ocessed widths fo r 
the Orion Nebula and W51 (see also iBell et HI 1201 If ). 
The left and right vertical lines show the 3a and la con- 
fidence limits from Figure [3 

4. DISCUSSION 

The multiple frequency switching technique reduces 
amplitudes of broad baseline variations and there- 
fore can be used to systematical l y "flatten" baselines 
(|Bell fc FeldmanI IT99TI IBeh et all [1992 IT991 . Com- 
pared to standard spline/polynomial baseline fitting, 
which allows for subjective decisions, the FS technique 
reduces baseline variations objectively. For spectral fea- 
tures that are much wider than the offset, h, FS can be 
likened to differentiation, which acts as a high-pass fil- 
ter in the frequency domain (see, e.g., [Hammin 
[Owen 1995). In this case, six overlaps approximates the 
sixth derivative with respect to frequency, which makes 
FS an efficient high-pass filter. This not only reduces 
broad variations in the spectrum, but in BASV's case 
(six overlaps, h < w ), it also reduces broad spectral line 
wings larger than the FS-offset, which manifests as line 
narrowing (s ee, e.g.. Figure [2l). 

We use the lLockman fc Brownlll975l model of the Orion 
Nebula and simulate the conditions of BASV's observa- 
tions at 6 GHz: telescope beam size, frequency range, 
channel width, offset, number of overlaps, and noise tem- 
perature. When modeling, we can avoid some difficulties 
that BASV faced when they processed the real spectrum. 
The real spectrum is abundant with spectral lines; after 
overlapping six times, the resulting spectrum is crowded 
with lines and reference lines. This situation is fur- 
ther complicated when significant noise fluctuations are 
present, which is BASV's case when An > 11, so clean- 
ing becomes a non-trivial and subjective procedure. In 
our case cleaning is trivial as we can model and clean 



each RRL individually, with and without noise. 

For all HnAn lines studied in BASV, the unprocessed 
width (i.e., before FS) is greater than the FS-offset used. 
We show that the processed width weakly depends on n 
and An (see Figure 21 the six FS-overlaps case). In fact, 
it remains about 0.8 /i for a wide range of n and An. 
This property of the FS technique makes it insensitive 
to changes in RRL widths when h < w. (It is for this 
reason that we justify the inclusion of BASV's processed 
widths for W51 at 6 GHz in Figure [H) The FS technique 
can be used to systematically eliminate broad baseline 
variations to detect weak lines, however, in BASV's case 
(six overlaps, h < w), this technique is not well suited 
to study RRL broadening. If one's purpose is to test the 
theory of spectral line broadening, the choice of the offset 
made in BASV {h < w) is not optimal. 

It is shown in the previous sections that FS reduces 
the S/N (see Equation ([6|)) and spectral line amplitudes 
(Figure [3]). As a result of amplitude reduction, the pro- 
cessed line temperature drops rapidly to the noise tem- 
perature Trms = 1 niK as n and An grow (the lower 
solid curve in Figure [5]). The model processed line am- 
plitude reaches 3 mK {3a level) at n = 224 (An — 11) 
and 1 mK {la level) at n = 241 (An = 14). The in- 
terval between n = 224 and 241 (between An = 11 and 
14) is a transition zone between observable spectral lines 
and noise. Figure |6] shows that the processed width be- 
havior in the transition zone imitates "narrowing" of the 
processed spectral lines observed by BASV in spectra of 
the Orion Nebula and W51 at 6 GHz. This "narrowing" 
is the result of "forcing" a program to fit Gaussians to 
features dominated by noise fluctuations — the procedure 
used by BASV. 

Spe ctra presented in iBelllllQQTl BASV and IBell et all 
120111 result from the overlapping and cleaning steps ex- 
plained in Section 3. Though BASV applied overlapping 
to the entire spectrum, cleaning was applied selectively 
to features located near the expected spectral line fre- 
quencies, even if the processed line temperatures were 
below 3 a. This subjective and selective approach could 
create "spectral lines" from random fluctuations at the 
< 3 cr level, which are abundant in the spectrum. 

For An > 14 the width errors predicted by our Monte 
Carlo simulation are comparable to the widths them- 
selves. Therefore, BASV's data points for An > 14 
shown in Figures [5] and [6] are likely to be misinterpreted 
noise features. 

In conclusion, FS is a useful technique for detecting 
and measuring weak spectral features, if required cor- 
rections are minimal, that is if the FS offset is greater 
than the line width: h > w. Based on our simulation re- 
sults, we argue that BASV's Trms ~ 1 mK observations 
of hydrogen RRLs from the Orion Nebula at 6 GHz are 
limited to An < 14. Within this limit, we find good 
agreement between their results and a Monte Carlo sim- 
ulation based on conventional Stark broadening theory 
and, therefore, we argue that BASV's findings do not 
necessitate a revision of RRL Stark broadening theory. 
We suggest that further tests of line broadening theory 
above this limit will require observations with sub-mK 
sensitivity and improved baseline stability. 

We thank Morley Bell for helpful communications 
when learning about FS. We thank Johannes Buchner 
for enlightening discussions about FS and various soft- 
ware tools. We thank Ron Maddalena for assistance with 
NRAO's 140 ft data archive. Finally, we thank Miller 
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